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This report is submitted to request the renewal of NASA Grant NAGW 1^0^ 
for the second year of funding, September 1, 1989-Oetober 31, 1990. 

I . Progress Report 

Our current efforts have focused on the dissociative recombination (DR) 
of Oj, a process of great importance in planetary atmospheres. This process 
is difficult to study experimentally because of the need to determine the 
dependence of the product electronic states and kinetic energies upon the 
vibrational distribution of the ion and the electron temperature. The 
knowledge of these characteristics of DR is needed to accurately model 
planetary ionospheres. Using a theoretical quantum chemical approach, we 
have studied in detail the generation of 0(*S) from DR, 

0+ + e" -> 0(*S) + 0( r D) , (1 ) 

where e - is an electron, both 0 atom products are excited, and the Oj 
molecular ion may be vibrationally excited. Building upon earlier 
studies 1 > 2 > 3 which showed that (1) occurs along a molecular state, we 
have now included the indirect dissociative recombination mechanism. The 
latter process is similar to (1) except that instead of direct electron 
capture into a repulsive molecular state leading to dissociation to atoms, 
electron capture first occurs into an intermediate vibrationally excited 
Rydberg state followed by predissociation into atoms along the repulsive 
state used in the direct mechanism, 


2 


0£ + e" -> 0 2 (n, v) ■> 0(*S) + 0(*D), (2) 

where 0 2 (n,v) is a neutral Rydberg state with principal quantum number n and 
vibrational level, v. These vibrationally excited Rydberg levels lie above 
the v=0 level of 0 2 and are resonances in the electron-ion continuum. When 
the energy of the incoming electron coincides with the position of the 
Rydberg vibrational level a marked perturbation occurs in the DR cross 
section. We have found that the DR rate can also be altered by these Rydberg 
states. 

The calculation of the full DR cross section including both the direct 
and indirect processes has required the calculation of the quantum defect of 
the 1 Zu Rydberg series and the autoionization width, both as a function of 
the internuclear distance. The cross sections have been calculated using 
molecular quantum defect theory. 4 Further details on these calculations are 
given in the first semi-annual progress report. 

Figure 1 shows the cross sections for both the direct and indirect DR 
processes. The lower curve at low electron energies is relatively smooth 
across the range of electron energies plotted and corresponds to the direct 
DR cross section, i.e. there are no Rydberg excited vibrational levels 
included in the calculation. However, the direct DR curve has several abrupt 
decreases at 0.23eV, 0.46eV, 0.68eV, and 0.90eV. These energies correspond 
to the positions of the v=1,2,3 and 4 vibrational ion levels above the v=0 
ion level. At each of these vibrational levels, a new channel for emitting 
the electron opens up and the DR cross section consequently drops. Some of 
these sudden decreases are partially obscured by the infinite series of 
resonances below each ion vibrational threshold which contribute to the 
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indirect DR process which is also plotted in the Figure. When the electron 
is emitted, the ion can be left in an excited vibrational level. The net 
result is electron impact vibrational excitation. Cross sections for the 
latter process have also been calculated and are discussed below. 

The full cross section curve describing both the direct and indirect 
process is also shown in Figure 1 . This curve has a very rich resonance 
structure giving rise to rapid variations in the cross section over narrow 
energy intervals. The full cross section curve is above the direct cross 
section curve even at low energies, where it is far from the first 
discernable resonance at 0.36eV. The source of the higher full cross section 
curve is the Rydberg level with principal quantum number n=6 and v=2. This 
level actually falls 0.01 47eV below the v=0 level of the ion. However it is 
broadened by predissociation by the dissociating state. The broadening 
allows a wing of this resonance to appear above threshold. This resonance 
plays a large role in increasing the DR cross section at typical ionospheric 
temperatures (see below). The first resonance appears at 0.036eV and 
corresponds to the n=9 Rydberg *1^ state in the v=1 vibrational level. The 
resonance has a typical Beutler-Fano profile corresponding to a positive 
Fano q, i.e. the cross section first decreases and then increases as the 
energy is increased through the resonance position. The next resonance for 
n=10,v=1 appears at 0.076eV with a similar shape to that for n=9. As the 
electron energy approaches the v=1 ion level at 0.23eV, the infinite series 
of v= 1 resonances interferes with resonances due to n=7,v=2 and n=4,v=6 at 
0.1 2eV leading to a complex cross section shape in this region. 

Figure 2 shows the DR rate constant for electron temperatures up to 
3000K. The upper curve is from the full DR calculation and the lower curve 





is for only direct DR in which the resonance states have been omitted. The 
full DR rate is 3.0 times larger that the direct DR rate at 300K and 2.6 
times greater near 800K. The increase is primarily due to the n=6, v=2 
resonance discussed above. 

Figure 3 shows the calculated vibrational excitation cross section 
corresponding to the following process: 

C> 2 (v= 0 ) + e“ -»• 0£(v-1) + e~. (3) 

The cross section for process (3) is calculated simultaneously with the 
calculation of the DR cross sections. Direct vibrational excitation by 
electron impact is often thought of as a slow process. However, its rate can 
be substantially increased if electronic resonance states can act as 
intermediate states. This is well known to be the case in the analogous 
excitation of neutral molecules where negative ion states play the 
intermediate role. There is very little in the literature on the excitation 
of molecular ions using neutral resonance states to assist the process. In 
the case of 0*, the cross section for excitation of the v=1 ion level has a 
sudden onset at the energy of the v=1 level, 0.23eV. Both the direct cross 
section, shown as the lower curve in Figure 3> and the full cross section, 
including both the direct and indirect processes, are shown. The direct 
cross section has abrupt breaks at the onset of higher excited vibrational 
levels. At the breaks, the cross section for v=1 excitation drops abruptly 
since a higher vibrational level can also be excited at these energies. At 
energies just below threshold, the cross section for the full DR process is 
considerably higher than that for direct DR and shows structure due to many 



5 


interfering resonance states. 

The vibrational excitation rates are shown in Figure H. The lower curve 
does not include the Rydberg resonances while the upper curve includes both 
direct and indirect DR. Clearly, the resonances play an important role in 
determining the magnitude of the excitation rate. The rate has a value of 
about 1.7 x 10 -11 cm3/sec at 800K. While this rate is not very large, it is 
important to remember that the repulsive state, which is studied here, 
does not cross the ion within the turning points of the v=1 vibrational 
level. However, there are other states which are known 2 * 3 * 5 to cross the ion 
within these turning points, e.g. 3 Iy and *A U . These states could lead to 
high vibrational excitation for v=1 . Studies of role of these states in 
vibrational excitation are currently in progress. 

II. Plan of Work for 9/1/89-10/31/90 

The plan of work for the next funding period will proceed in accord with 
the research discussed in the original proposal. The plan includes the 
completion of the ab initio calculation of the Rydberg and valence potential 
curves needed for the study of both direct and indirect DR leading to 3 P 
atoms. Only states having significant electron capture widths will be 
determined. These states will be used to calculate cross sections and rates 
for DR leading to 3 P atoms, using techniques developed during the first 
year of funding. The electronic widths of states leading to *D atoms will be 
calculated. For those states with large widths, the potential curves will 
also be determined. The important Rydberg states contributing to indirect 
recombination in the generation of *D atoms will be determined. DR cross 
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sections and rates for generating *D atoms will be calculated for both the 
direct and indirect processes. The results will yield a detailed picture of 
the DR of Oj and will be invaluable for the modelling of the chemistry of 
planetary ionospheres. The kinetic energies and quantum yields of excited 
and ground state atoms will be determined as a function of electron 
temperature and ion vibrational excitation. 
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AB INITIO STUDIES OF DISSOCIATIVE RECOMBINATION 
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ABSTRACT 

Quantum chemical calculations of the dissociative 
recombination of Oj and N^ are reported. An approach for 
calculating auto ionization widths from high principal 
quantum number Rydberg states is summarized and an 
example is presented for the lowest s Hg dissociative 
state of 0 2 . For Oj, the 'Ey state is the sole source of 
O^S) from the lowest 10 vibrational levels of the ion. 
For OC'D), the B 3 !^, *1^, and 1 *A U states are the sources 
of 0( 1 D) from the lowest three vibrational levels. Rate 
coefficients for generating 0(*S) and O^D) at ionospher- 
ic temperatures are reported. For Nj, states of 3 n u , *Eg 
and 3 Ag symmetries intersect the ion at the v-0 level. 
The configurational structure of these states indicates 
that the 3 n u states are likely to be the most important 
routes for dissociative recombination from vO. 


1 . INTRODUCTION 

Interest in dissociative recombination (DR) first arose out of the 
need to understand the properties of the Earth's ionosphere. 1 ^ 
Initial attention focused on the DR of O 2 * Nj and N0 + . While there has 
been much progress in the measurement of the total DR rate coefficients 
of these ions, the dependence of the quantum yield of excited and 
ground state atoms upon the electron temperature and ion vibrational 
level is only currently being unraveled. Impediments to experimental 
progress have included the difficulty of generating ions in specific 
vibrational levels and simultaneously identifying the atomic products. 
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Nevertheless there has been recent progress toward this goal. 3] »^] *5] 
An alternative theoretical approach based on spectroscopic identifica- 
tion of the potential curves which provide routes for DR is difficult 
since these states often do not have dipole allowed transitions from 
the ground state. Furthermore, electron capture occurs at high energies 
on the repulsive wall, well above the dissociation asymptote and out of 
range of the RKR approach. Additional difficulties include the paucity 
of information on electron capture widths although progress has been 
made in spectroscopic analysis of R y d b e r g - v a 1 e n c e 
interactions^ »7] >8] .9] from which widths can be derived. 

Here, I review recent theoretical ab initio results obtained in this 
laboratory for the DR of 0 Z and describe some preliminary results on 
the DR of N z . For 0 Z , potential curves, electron capture widths and 
rate constants have been obtained for generating *S and 'D atoms, the 
upper states of the green and red lines respectively. For N z , poten- 
tial curves for DR of the v=0 ion level have been computed. 

2. WAVE FUNCTIONS 

The wave functions calculated here are based on orbitals determined 
in either multiconfiguration or complete active space self consistent 
field (MCSCF or CASSCF) calculations. Each of the orbitals, <J> j , is 
expanded in a large set of Gaussian basis functions, Xi» 

♦j * l c j , i Xi (1) 

Some of the Xi are individual Gaussian primitives 1 ^] while others are 
fixed combinations (contractions) of Gaussian basis functions deter- 
mined from atomic calculations. 11 ] The orbitals are used to construct 
large scale configuration interaction (Cl) wave functions containing 
many configurations, each of which describes a possible distribution of 
the electrons in the molecular orbitals for a particular state symme- 
try. For example, for the 1 Eg ground state of N z , the most important 
configuration near the equilibrium separation is: 



( 1 a g ) 2 (1a u ) 2 (20g) 2 (2o u ) 2 (3a g ) 2 (Wu)‘' 


( 2 ) 


The la and 2 a orbitals are mostly N atom Is and 2s orbitals, respec- 
tively. The 3tJg orbital is a bonding combination consisting mostly of 
the N atom 2p orbitals lying along the bond axis. The 1ir u orbitals are 
out of plane bond orbitals consisting mostly of the N 2p orbitals that 
are perpendicular to the bond axis. Each configuration can have one or 
more spin couplings (referred to below as terms) for a given total spin 
symmetry. A wave function is constructed consisting of a superposition 
of many terms, $j[, each of which has the same total spin and spatial 
symmetry, 

N 

*j = l Cj.i *i- (3) 
i=1 

In the CASSCF^] approach a set of active orbitals is defined and 
N(CASSCF) terms are generated in (3) by taking all possible excitations 
of the electrons within this set while restricting each term to have 
the total spin and spatial symmetry of interest. Both the cj^ and the 
C j t i are determined by variationally optimizing the total energy of j . 
The MCSCF approach is identical except that usually a subset of the 
full active space is used. Since the orbitals are expanded in large 
basis sets, Xi* the number of orbitals used in the CASSCF is only a 
subset of the orbitals that can be formed from the basis set. The Cl 
wave functions are constructed by exciting electrons out of a small set 
of reference configurations into the remaining (virtual) orbitals. The 
reference configurations often consist of all the configurations needed 
to properly dissociate that particular state to the proper separated 
atoms limit. This approach leads to wave functions that are not biased 
to favor particular internuclear distances. Some of the wave functions 
reported here for N 2 include the 2o orbitals in the active space and 
use the entire CASSCF wave function as a reference set. Taking at most 
one or two electrons from each of the reference configurations and 
promoting them to the virtual and other valence orbitals leads to a 
singles and doubles Cl denoted CISD. The Cl wave functions are expanded 
over the CASSCF orbitals and have the form of (3) with 


N(CI)>>N(CASSCF) . Some of calculations used here employ large scale 
wave functions having more than 100,000 terms with the Cl coefficients 
determined by direct Cl techniques. 1 33 Calculations have recently been 
completed on several states of N 2 and 0 2 which provide examples of the 
sort of accuracy which can be expected with these techniques. For the 
ground state of N 2 , 14 reference configurations (needed to properly 
dissociate the ground state to l *S atoms) are used in the Cl with a 
basis set consisting of 6s,3p,2d and If type contracted Gaussians 
leading to a Cl wave function having 176,536 terms. The calculated 
(experimental 1 ^) spectroscopic constants are 2340. 57(2356. 56cm -1 ) for 
(jo 0 , 1 4.27(1 4.3244cm -1 ) for u e x e , and 2. 0856(2. 074347a 0 ) for the 
equilibrium distance, R e . Similar calculations on the ground state of 
§2 using a basis set of the same size leads to a 228,036 term wave 
function. The calculated (experimental 1 5]) results are 1565.47(1580.19) 
for u) e , 10.93(1 1 .98cm -1 ) for w e x e , and 2. 2999(2. 28l9a 0 ) for R e . For 0 2 , 
a series of vertical excitation energies have been calculated. The 
calculated (experimentally derived) results are 5. 8982(5. 90±0.05eV) for 
‘ZuT 6. 1318(6. 11±0.05eV) for 3 A U , 6. 2766(6. 19±0.05eV) for s Eu* and 
8.66(8.6l±0.05eV 1 ) for 3 E U 7 All the experimental energies are derived 
from Ref. 15 (except for 3 Zu by extending the RKR results with an ab 
initio repulsive wall. 

Recent results have been obtained for N z using a [43,3p,2d, 1 f] 
contracted Gaussian basis set and a CISD reference set consisting of 
the CASSCF configurations with 2a in the active space. For the C 3 n u 
state the calculated (experimental ) 1 ^ ^ results are 1 1 .01 (1 1 ,05eV) for 
T e , 1 986(2047cm -1 ) for u> e , 25 . 05 ( 28 . 44 cm -1 ) for aj e x e , and 
2 . 1 932 (2 . 1 707a o ) for R e . For the w*A u state the results are 
8.83(8.84eV) for T e , 1532(1559cm -1 ) for u) e , 1 1 .67( 1 2.01 cm -1 ) for u> e x e , 
and 2. 4192(2. 3977a 0 ) for R e . The highly accurate values of the excita- 
tion energies and the other spectroscopic constants give U3 confidence 
that these repulsive curves will be located quite accurately relative 
to the ion potential curve. 

3. POTENTIAL CURVES FOR 0 2 

Potential curves describing DR in 0 2 and based on a C3s,2p,1d] 



contracted Gaussian basis set and first order Cl wave functions have 
been described previously. 1 ^ These calculations showed that the ‘Zy 
state is the only state that can generate O^S) from the low vibration- 
al levels. It crosses the ion between the outer turning points of the 
v=1 and v=2 vibrational levels and dissociates to the J D + *S asymp- 
tote. The next accessible state is 5 3 n g which crosses the outer turning 
point of v=10. Among the states dissociating to OC'S) there are three 
routes which cross the a‘*n u metastable state between the turning points 
for v=0. These routes with their dissociation limits shown in parenthe- 
ses are. 5 3 lIg( 3 P + l S), 5 3 n u ( 3 P + ‘S), and 4 1 Tig C 1 D + *S). Because of its 
spin symmetry, the 4 1 n g state will have only a small DR rate coeffi- 
cient from the metastable state. The remaining two states can lead to 
hot 0 atoms, each with 3.4eV kinetic energy. Future experiments 
designed to measure the l S quantum yield must be careful to be certain 
that no a‘*n u 0* is generated which would interfere with the determina- 
tion of the *S generated from the ion ground state. 

For generating CK'D), in addition to the 'Ey state, the B 3 Zu state 
crosses the ion near the inner turning point of v»0 and leads to 3 P + 
*D atoms. A l A u state leading to two l D atoms crosses the ion near the 
outer turning point of v=0. The lowest valence state of 1 Ilg symmetry 
crosses the ion near the inner turning point of v=0 and can lead to *D 
atoms via a curve crossing with the 2*1^ state. At higher energies, the 
lowest 3 H g state crosses the ion near the inner turning point of the 
v=*3 ion level and can lead to *D atoms via an avoided crossing with the 
2 3 Ug state. The next accessible states for *D production are 2 3 E g ( 1 D + 
3 P), 2 3 Ilg ( 1 D + 3 P), and 2 3 H U ( 1 D + 3 P) near the outer turning point of 
v= 3> 3 3 n g (*D + 3 P) and 2 1 n g ( l D + *D) near the outer turning point of 
v=4. As has been pointed out earlier 1 several of these routes are 
expected to have small electron capture widths. 

The reader is referred to Ref. 17 for a discussion of additional 
routes that lead to only ground state atoms. 

The important B 3 Zy, *A U , and ‘E^ states identified above have been 
the subject of larger scale calculations involving [6s,3p»2d,1 f ] 
contracted Gaussian basis sets and Cl wave functions that include all 
single and double excitations to the virtual space. 1 > ^3 Second order 



Fig. i. Dissociative first order (dashed) and second order 
(solid) potentials for 0 2 plotted with the ion potential 
from Ref. 15. 


curves for the Hg states have not been obtained because of their small 
widths (see below). Both the earlier curves and the large scale curves 
are plotted in Figure 1. The ion curve is the RKR^5] curve plotted at 
the experimental ionization potential above the large scale calculated 
ground state. The ion curve has been shifted to larger R to compensate 
for the difference between the calculated and experimental ground state 
R e . In order to plot the first order curves on the same figure they 
have been shifted to larger R to correct for the difference between the 
fi?st order and second order calculated R e 's. Also, the first order 
curves are shifted in energy so as to be plotted relative to the first 
order ground state. Relative to the ion, after the shift in the energy 
and distance coordinates, the curves are qualitatively similar. The 
largest difference is for *A U for which the second order curve is 
shifted by about 0.02a o to smaller R near v=0 relative to the first 
order rusult. 


*J. ELECTRON CAPTURE WIDTHS FOR 0 2 

The DR cross section is approximately directly proportional to the 


electron capture or autoionization width, r, a matrix element of the 
Hamiltonian operator, H, given by Fermi's golden rule: 

T(R) = 2irp | <'Ff | H | | 2 . (4) 

Tj is the wave function of the properly antisymmetrized product of a 
"free" electron in a coulomb like orbital and the wave function for the 
ion. Tf is the electronic wave function for the neutral state which 
describes the DR products and p is the density of states. Rydberg 
orbitals with high principal quantum numbers, n, have been used to 
describe the coulomb like orbital. Except for the Rydberg orbital, all 
the orbitals in (4) are bound and have large amplitude only near the 
molecule. As a result, it is only important to know the amplitude of 
the Rydberg orbital near the molecule. The inner part of coulomb 
orbitals are well known to be very similar to the inner part of high 
Rydberg orbitals. 20 ^ With the coulomb orbital represented as a high 
Rydberg orbital all the orbitals in (4) are bound and are represented 
in terms of Gaussian basis functions. 

For n*>»l, where 8. is the angular momentum quantum number, the 
matrix element in (4) can be closely represented by 

| <4'f | H | 'i'5 > = ( 1 /n*3)kexp(-c/n#) (5) 

where k and c are constants determined from the calculation of the 
matrix elements on the left side of (5) for the highest Rydberg states. 
These Rydberg states are nearly hydrogenic and the energies can be 
represented as E=-(1 /2(n*) 2 ) where n*=n-S and 6 is the quantum defect. 
The density of states can be written as 

p - 1/(E(n*-1/2) - E(n*+1 /2) ) (6) 

which can be closely represented as 

p = n*3exp(-1 /(2n* 2 ) ) . (7) 


Substituting (7) and (5) into (4) gives 



r - 2Trkexp(-(c/n*)-( . 5/n* 2 ) ) . 


( 8 ) 


Taking the limit of n* + ® gives the threshold capture width, 

r = 2irk. (9) 


The Rydberg wave functions have been refined further by optimizing them 
in the field of an optical potential due to the states. Feshbach 
projection operators 21 ^ have been defined which partition the total 
space such that V - PY + QV where PK contains the Rydberg states to be 
optimized and Q¥ are the Yf states described above. In the usual 
Feshbach projection operator formalism, an eigenvalue equation is 
derived for PV that is difficult to solve because it contains an 
energy-dependent optical potential. However, using partitioning 
techniques, 22 -! P 1 ? can easily be determined by solving the usual Cl 
problem. Writing the total wave function as ¥ = P¥ + QV, we can write 
the Schroedinger equation in matrix form as 



( 10 ) 


where Hpp = PHP, Hqp = QHP, etc. 

Multiplying the matrices in (10) leads to: 


HppP'i' + HpqQH' = EPH 1 

and 

HQpPf + HqqQ'P = EQV. 
From Eq. (12) we have, 


(ID 

( 12 ) 


Q¥ - HqpP'i'/ (E - Hqq). 


(13) 


Substituting Eq.(13) into Eq.(11) leads to a matrix optical potential 
for P¥, 


(H PP + H P qHq P /(E - Hqq)) Pf - EP¥. 


( 1*0 



It is difficult to solve (14) directly for PT since E is on both sides 
of the equation. Nevertheless, PY can be determined by retaining the 
coefficients of the P^ terms resulting from the diagonalization of the 
H matrix in Eq.(10). Before diagonalizing H we must project out the low 
energy QT roots in order to prevent them from mixing into P'C. There- 
fore, it is necessary to first solve for the O'? roots by diagonalizing 
H qq . The Hqq portion of the H matrix in Eq. (10) is transformed to 
project out the physically meaningful Q'F roots. P*P is then determined 
by diagonalizing the new H matrix with the transformed Hqq. The net 
effect is to provide additional correlation to the P^ space from the 
nonphysical roots in the Q space. P^’s for successive principal quantum 
numbers are used to represent in Eq.(4). The widths obtained by this 
procedure are then extrapolated to the continuum as described above. 
The importance of the correlation added to the P space from the higher 
terms in Q space has been emphasized by Hazi who has described an 
approach similar to that used here. 2 ^^ 

An application of this approach to the calculation of the widths of 
the B and L 2 II states of NO has already been discussed. 2 ^ An 
additional example of this approach is given in Table 1 . for the 0 2 3 JI g 
width. This is the only 0 2 width for which there are experimentally 
derived matrix elements for comparison. An analysis of the released 
kinetic energy line width for the v=1 level of the C 3 IIg state of 0 2 
indicates that the electronic matrix element between the C Rydberg 
state and the lowest valence 3 n g state is 0.079eV.®^ An additional 
determination based on the experimental line widths of the lowest 4 
vibrational levels of the C state obtains a matrix element of 
0. 0625(+0. 006,-0. 004 )eV. 9 ^ 

All calculations in Table 1. have been done at R=2.28l9a 0 . A 
[3s,2p,1d] basis set on each center was supplemented with a set a 18 
diffuse s Gaussians at the midpoint and the Rydberg orbitals were 
determined by the Improved Virtual Orbital method. 2 ^ The 3 Hg valence 
state was described with a full valence Cl consisting of 46 terms. For 
the C state, a valence Cl on the ion ground state was coupled with the 
Rydberg orbital and each of the three valence virtual a g orbitals 
leading to a 636 term Cl. Using the procedure outlined above leads to 


the widths shown in Table 1 . Using the widths for n-8,9 gives the 
expression, r- 0.00217exp(1 .39/n#~(0.5/n* 2 ) ) . The extrapolated width 
for n*-»“ is 0.00217eV. This width is about two orders of magnitude 
smaller than the widths for the three important dissociative routes 
identified above. As a result, the 3 n g state was not included in the DR 
rate coefficient calculations discussed below. 


Table 1. Electron-Ion Capture Widths (r) from 3 n g Rydberg States 


n 

n* 

p ( eV 1 ) 

r(eV) 

9 

7.8897 

17.903 

0.00257 

8 

6.8896 

11 .892 

0.00263 

7 

5.8902 

7.4022 

0.00272 

6 

4.8913 

4.2110 

0.00285 

5 

3.8930 

2.0972 

0.00311 

4 

2.8959 

0.84005 

0.00367 

3 

1 .8984 

0.21774 

0.00522 


The matrix element calculated here for n=3 is 0.06l7eV. Since the 
electronic matrix element varies with R, the experimentally derived 
matrix element is an averaged value. However, since the variation of 
the matrix element with R over the relevant range of nuclear vibrations 
is estimated to be less than about 10$, the results show that rela- 
tively small wave functions can yield interaction matrix elements that 
are in good agreement with experiment. 

5. DR RATE COEFFICIENTS FOR 0 2 

The expression for the direct DR cross section derived by Giusti 2 ^ 
has been used to calculate cross sections and rate coefficients for 
generation of OC'S) and 0( 1 D) along the B 3 £^, and 1 L U diabatic 
potential curves described above. For DR leading to 0( X S ) , "* ® the v»2 
level of the ion has the highest rate coefficient for the lowest eight 
vibrational levels and is a factor of 78 greater than the coefficient 
from v=0 at an electron temperature (T e ) of 300K. A comparison of the 
calculated rate coefficients to those deduced from atmospheric models 
shows that the atmospheric models apply to molecular ions which are not 



in a Boltzmann distribution but which are instead significantly 
populated in the v>0 levels. In contrast to the situation for 0(*S) 
generation, the total rate coefficient 1 ^ for producing 0(*D) is 
greatest from v=0 and is a factor of 1 .9 greater than the coefficient 
from v=2 at T e =300K. The rate coefficient for v=»0 i3 in good agreement 
with those deduced from atmospheric models. The calculations show that 
the dominant route for producing 0(*D) from v=0 is *A U and 3 Zy is the 


Table 2. DR Rate Coefficients for Production of 0( l S) for 600£T e S1000K. 


Ion Vibrational Level 

Rate Coefficient(cm 3 /sec) a 


0 

0.24(+.12,-.06)[-0.14] 


1 

5.0( + 1 .2,-1 .DH-0.36] 


2 

15.(+2.,-3.)[-0.50] 


3 

*1.8( + 3.3,-2.3)[-0.77] 


4 

4 .9 (+0 .6 , -1 .2)[-0.32] 


5 

1 .2(+1 .9, -0.9) [-0.90] 


6 

5.1 (+0.6.-1 .1 )[-0. 56] 


7 

1 .4(+0.9,-0.9)[+0.06] 



a. All rate coefficients have been multiplied by 10®. The exponent of 
the temperature dependence is in square brackets. The rate coefficient 
for v=0 is 2 .4( + 1 .2 ,-0.8 )x1 0 -1 °x( (T e /800)~° • 14 )cm 3 /sec. 


dominant route from v=1 and v=2. 

Since the electron temperatures in ionospheric models are near 800K, 
we list in Table 2. the rate coefficients for producing CK'S) for 
600ST e £1000K for each of the lowest eight ion vibrational levels. 


Table 3. DR Rate Coefficients for Production of 0(*D) for 600£T e £1000K. 


Ion Vibrational Level 

Rate Coeff icient(cm 3 /sec) a 

0 

1 .37(+.l4,-.l6)[-0.51] 

1 

1 .09 (+.12, -.14) [-0.55] 

2 

0.78(+.13,-.l6)[-0.4l] 


a. All rate coefficients have been multiplied by 10 7 . The exponent of 
the temperature dependence is in square brackets. The rate coefficient 
for v=0 is 1 .37( + .1 4 ,-. 1 6)x1 0 -7 x( (T e /800)~ o • 5 1 )cm 3 /sec. 


Because of the position of the intersection of the dissociative curve 
with the v=5 level, the coefficient from v=5 is highly sensitive to an 
energetic displacement of the dissociative route within its expected 
uncertainty (O.leV). Table 3 has the 0( 1 D) rate coefficients in the 


same temperature range from each of the lowest three ion vibrational 
levels. 

6. DISSOCIATIVE ROUTES FOR N 2 

Preliminary results are reported here which identify the important 
dissociative routes for the v=0 level of the ground state of N 2 . These 
results are based on large scale second order Cl calculations using 
atomic natural orbital basis sets. 27] The details will be reported 
separately. 

The dominant dissociative routes from v=0, 2 1 I g ( 2 D + 2 D), C ,, n u (‘*S + 
2 P). 3 3 rt u ( 2 D + 2 D), and G 3 A g C , S + 2 D), are shown in Fig. 2. 

The dominant configuration of the state is 

. . . . ( 3c g ) 2 ( 1 iTux ^ ( l^uy) Oifgx) ^ ^ ^ gy ^ 

where .... denotes ( 1 o g ) 2 ( 1 a u ) 2 (2o g ) 2 (2a u ) 2 . This configuration differs 
by a triple excitation from the dominant configuration of the ion plus 
a "free" or Rydberg electron, 

....(3o g ) (1tr ux ) 2 (1ir U y) 2 (RYDo g ). 

The matrix element of the Hamiltonian between these two configurations 
is zero and these dominant terms will not contribute to the electron 
capture width. The magnitude of the width will be determined by matrix 
elements between secondary terms in both states and may be small. As a 
result, 2 2 E g may make only a small contribution to the total DR rate 
coefficient from v=0. 

The G state also passes through the v=0 level of the ion and 
dissociates to H S + 2 D giving atoms with 1 .72eV kinetic energy from 
v=0. The dominant configurations of this state, 

....(3c g ) 2 (lTTy X ) ( l^uy) ( ^ gx^ f 
• . . . (3 o g ) 2 ( It^ux ^ ( l^uy) ( 1lT gy) > 

• •.•(3Cg) 2 (lTTy X ) ( ^ gX ^ ( iTTgy) > 

• • ■ • (30g) 2 ( 1 1T gx ) ( 1 TT gx ) ( ^ TT gy ) , 



differ from the ion plus a continuum electron by a triple excitation. 
Therefore, as for 2 z Eg, the dominant configurations will not contribute 



INTERNUCLEAR DISTANCE (Bohr) 

Fig. 2 Dissociative routes for the v=0 level of Nj. 

The ion potential is from Ref. Ht. 

to the width and this state is likely to have a small width and a small 
DR rate coefficient. This is an important conclusion for the escape of 
N atoms from Mars since this state and possibly the C* state (through 
an avoided crossing with the C state) are the only states that can 
supply enough kinetic energy for escape from v=0. 

The primary configuration of the C' and 3 3 Hu states, 

. . . . (30g) Oifux^d^uy) ^^gx^ (I'fgy) > 

differs by a double excitation from the ion ground state plus a 
continuum electron. The capture width for this state may be larger than 
that for 2*Eg and G 3 A g and it appears that these 3 II U states may be the 
dominant states for DR from v=0. The importance of the 3 n u states is in 


good agreement with the results of Michels. 28 ^ The C' state has an 
avoided crossing with the lower C s n u state which dissociates to 1, S + 2 D 
atoms. The determination of the states of the product atoms due to 
initial dissociation along the C' state is currently under study. 

7. THE ROLE OF RYDBERG STATES IN DIRECT DR OF N 2 

Rydberg state potential curves having the ground state of the ion as 
core cannot cross the ground state of the ion and in addition have very 
small electron capture widths. These states are usually not important 
for direct DR electron capture. However, because the first excited 
state of N 2 is only 1.1 eV above the ground state, Rydberg states having 
the excited state core can cross the ground state of the ion and 
provide routes for DR. 

For 1 Eg symmetry, an N 2 Rydberg state can be formed by binding a ir u 
Rydberg orbital to the excited A 2 n u state of N 2 . This state will have 
an avoided crossing with the 2‘Eg state described above and can provide 
a channel for DR. A comparison to analogous states in 0 2 indicates 
that the 4pir ‘Eg state having the A core will pass through the ground 
state of N 2 near the equilibrium separation. If the 4pir ‘Eg state 
should prove to have a large capture width, this route could be an 
important route for DR in N 2 . The calculation of the potential curve 
for this state and the electron capture width are in progress. While 
the Rydberg states are diffuse, their primary configurations differ 
from the ion by a double excitation and they may have larger electronic 
widths than the ‘Eg valence state which has primary terms that differ 
from the ion plus a continuum electron by a triple excitation. It 
appears that a multistate treatment involving both the valence states 
and the excited core Rydberg states may be needed to calculate the 
cross sections for DR of N 2 . 
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